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Abstract. For systems in an externally controllable time-dependent potential, the 
optimal protocol minimizes the mean work spent in a finite-time transition between 
given initial and final values of a control parameter. For an initially thermalizcd 
ensemble, we consider both Hamiltonian evolution for classical systems and Schrodinger 
evolution for quantum systems. In both cases, we show that for harmonic potentials, 
the optimal work is given by the adiabatic work even in the limit of short transition 
times. This result is counter-intuitive because the adiabatic work is substantially 
smaller than the work for an instantaneous jump. We also perform numerical 
calculations of the optimal protocol for Hamiltonian dynamics in an anharmonic 
quartic potential. For a two-level spin system, we give examples where the adiabatic 
work can be reached in either a finite or an arbitrarily short transition time depending 
on the allowed parameter space. 
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1. Introduction 

Recently, there has been considerable progress in the thermodynamic description of 
small (bio-) systems [U [2] which are prone to thermal fluctuations and typically driven 
far out of equilibrium. It has been shown that thermodynamic quantities obey various 
exact work and fluctuation theorems [31 HI El El Ej. The Jarzynski relation connects 
the equilibrium free energy difference with an average of work values obtained from a 
series of nonequilibrium transitions [51 E]. While biomolecules are typically described 
by a Langevin equation where effects of the thermal bath are included in friction and 
fluctuation terms, the Jarzynski relation has originally been derived in the framework 
of purely Hamiltonian dynamics with initial canonical equilibrium distribution [5]. 
This situation corresponds to an experiment where the heat bath of temperature T 
is decoupled from a small Hamiltonian system at time r = 0. During a time intevall 
r G [0, t] a control parameter A is varied time-dependently. Such dynamics has been 
used in various generalizations of the Jarzynski relation pi [TO] and for the optimization 
of the free energy reconstruction from nonequilibrium work data [TTI [12| [T3j . 

While this approach is elegant from a theoretical point of view, it is important 
to note that such a dynamics has quite different properties than a system subject to 
a permanent heat bath modelled e. g. by a Langevin equation. The most impor- 
tant difference concerns the quasistatic work which is achieved for infinitesimally slow 
transitions t — > oo, independent of the detailed shape of the protocol for the control 
parameter A. For a system coupled permanently to a heat bath conserving the canonical 
distribution for fixed control parameter, this quasistatic work is equivalent to the free 
energy difference. This property gives rise to the free energy calculation method termed 
thermodynamic integration where the free energy is approximated by the work for a 
very slow transition. For purely Hamiltonian dynamics, however, the quasistatic work 
(which then is called adiabatic work) can exceed the free energy difference p2]. Still, 
the Jarzynski relation is valid and can be used to estimate free energy differences. It 
can be shown that under reasonable assumptions (mainly ergodicity), both the work is 
a monotonically decreasing function of the transition time t for a given protocol shape 
and the adiabatic work is independent of the shape of the protocol [H]. 

There has been considerable interest in deriving fluctuation theorems and analogues 
of the Jarzynski relation also for systems in the quantum domain [151 EH EL71 (HI [El [20]. 
Here, it is more subtle than in the classical case to define thermodynamic quantities 
such as work on a single trajectory [TTJ, [181 US]- Apart from recent progresses for open 
quantum systems [19] , most studies [T5l ED, [22], [23] consider the case where the energy 
is measured at beginning and end, defining the work as the difference between final and 
initial energy. These studies rely on Schrddinger dynamics which has similar properties 
as Hamiltonian dynamics. In particular, the quasistatic work then is also different from 
the free energy difference. 

In this paper, we study optimal protocols of the control parameter A(r) leading to a 
minimal mean work for a finite given transition time t. This question has been tackled 
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for Langevin dynamics in harmonic potentials [2H [25] with the result of unexpected 
singularities such as jumps in the overdamped case and even delta-type singularities in 
the underdamped case. For purely Hamiltonian or Schrddinger dynamics in harmonic 
potentials, we will show that the optimal protocol is degenerate and that the optimal 
work is equal to the adiabatic work. This result implies that the adiabatic work, which 
is achieved for an infinitely long transition time for any protocol, can be reached in an 
arbitrarily short time by using an appropriate protocol for the control parameter. This 
surprising result shows that purely Hamiltonian dynamics can beat Langevin evolution 
since for the latter short transition times, t — > 0, yield the work for an instantaneous 
jump. In the totally different context of collective escape over a barrier, it has been 
shown that Hamiltonian dynamics also beats Langevin dynamics in the sense that it 
leads to a faster escape [26]. We will also calculate the optimal work numerically for 
an anharmonic potential. However, it is difficult to decide from the numerical data, 
whether the adiabatic work can be achieved in a finite transition time also for such an 
anharmonic potential. 

2. Hamiltonian dynamics 

In this section, we will study minimum work protocols for Hamiltonian dynamics in 
three different types of time-dependent potentials: (i) a moving harmonic potential, 
(ii) a harmonic potential with time-dependent stiffness and (iii) an anharmonic quartic 
potential. 

2.1. Moving harmonic potential 

We consider a particle of mass m subject to a moving harmonic potential 

V(x,r) = ^(x-\(r)) 2 , (1) 

where k is the (constant) stiffness of the potential. The minimum of the potential 
A(t) is changed time-dependently from an initial position Aj = to a final position A/. 
The Hamiltonian dynamics of position and momentum of the particle at time r then is 
governed by the equations of motion 

x = p/m, 

p=-V'(x,T), (2) 

where the prime denotes the derivative with respect to x and the dot the derivative with 
respect to the time r. The particle initially is in thermal equilibrium with probability 
density 

p(x , Po ) =Afexp[-(3V(x ,0)]exp[-f3p 2 /(2m)] (3) 

for the initial position x$ and the initial momentum po with normalization constant M 
and inverse temperature (3 = 1/T with Boltzmann's constant ks = 1. 
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Since there is no heat transfer during the (purely Hamiltonian) transition, the work 
is given by the change of internal energy. Averaging over all initial conditions yields the 
mean work 



W 



2m 

(p!> 

2m 2 



(V) 
k 



((x 2 ) - 2A (x) + A 2 ) 



(4) 



where (•) denotes the average over the initial distribution of position and momentum 
(J3]). It can easily be shown that the variances (x 2 ) — (x) 2 and (p 2 ) — (p) 2 are time 
independent for any protocol A(r). Thus, the work (j4j) can be written as 

v2 , 1* 



2m 



+ - ((x) 2 - 2A (x) + A 2 ) 
z Jo 



which becomes zero, i.e., minimal for 

(P(*)> = o, 
(x(t)> =A/. 

The evolution equations for these mean values are given by 
d 



(5) 



(6) 



dr<*> 



(p) /m, 
- k((x) - A) 



(7) 



with initial values (x(0)) = and (p(0)) = 0. With the freedom to choose the continuous 
function A(r) which corresponds to an infinite number of free parameters, it is quite 
obvious that the two final conditions can easily be met. In fact, the optimal protocol is 
highly degenerate. 

For a specific choice, taking, e. g., a third order polynomial 



A(r) = a(r/t) + b{r/t) 2 + (A/ - a - b){r/tf 



(8) 



with two free parameters a, b suffices to meet the optimality conditions. The solution of 
the evolution equations (JTj) can be solved analytically and, inserted into equation (J6]), 
leads to 

12A/ sin {ut/2) - 6cut\ f cos (cut/2) 



(12 - ujH 2 ) sin (ut/2) - 6cut cos (cot/2) ' 

-3cu 2 t 2 Xf sin (cut/2) 
[12 - co 2 t 2 ) sin (cot/2) - Qcot cos (cot/2) 



(9) 



with co = Jk/m. For short transitions times iCl, the coefficients diverge as a ~ 1/t 2 
and b ~ -1/t 2 . 
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As a quite different choice, both conditions can also be met by a linear protocol 
with delta singularities at the boundaries, in analogy to the underdamped stochastic 
dynamics case [25]. Indeed, it can easily be verified that the protocol 

y(r) = X f j + ^[5(r)-5(r-t)] (10) 

fulfills the conditions ([6]). This expression could also have been obtained as the limit of 
the underdamped optimal protocol for a vanishing friction coefficient 7^0. 

Thus, for any given total transition time t the minimal mean work is given by 

W* = 0. (11) 

For infinitely long transition times, quite generally, the work is given by the adiabatic 
work. For the studied model system, we have W &d = 0. As shown above, this 
lower bound on the work can even be achieved for arbitrarily short transition time. 
This surprising result shows that purely Hamiltonian dynamics can beat Langevin 
evolution where short transition times t —>■ yield the work for an instantaneous jump 
W ip = kXj/2. 

The optimal protocols show a distinct positive peak in the first part of the transition 
and a distinct negative peak in the second part of the transition. This peak structure 
can be understood intuitively in complete analogy to the underdamped dynamics case 
[25] . The first peak accelerates the particle and leads to a finite mean velocity which 
is necessary to reach a final mean particle position at the potential minimum. The 
second negative peak decelerates the particle, thus recovering the invested work from 
the kinetic energy of the particle. 



2.2. Harmonic potential with time dependent stiffness 
We consider the time dependent potential 
V(x,r) 



Mr)_ x2 



(12) 



where the stiffness A(r) is changed from an initial value A(0) = Aj to a final value 
A(t) = A/ in a finite time t. The work exerted on the system during the finite time 
transition is given by the change in internal energy 



W 



2m 



+ (V) 



(P 2 ) , A 



2m 



— { x 
2 



J 



(13) 



The dynamics ([2]) leads to coupled evolution equations for the mean squared position 
w = (x 2 ) and the mean squared momentum z = (p 2 ) 

z = — mXw, (14) 

(15) 



w 



2z 2 , 
— Xw. 

m z m 



Multiplication of equation (fl5"l) with a factor w and insertion of equation (fill) leads to 
the relation 

ww = — - (zw + wz) (16) 
m 2 
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which, using the equilibrium initial conditions w(0) = T/\, w(0) = 0, and z(0) = mT 
can be integrated to yield 

m 2 2 mT 2 

zw = —-w H : — . (17) 

4 Ai 

At the final time t, minimization requires w(t) = which is equivalent to the final 
condition 

„ , , mT 2 , . 

z*(t) = - — . 18 

The (mean) total work thus becomes 
1 



W 2 



\ f w(t)-T+-^—-T 
w{t)Xi 



(19) 



Optimizing this expression with respect to w(t) leads to 

w*(t) = (20) 
\J\Xf 

which yields the optimal work 

W* = T\J^-1\. (21) 




Similar to the first case study, the optimal protocol is highly degenerate. For a third 
order polynomial, the equations of motion cannot be solved analytically. However, as 
expected, the two parameters can be chosen such that the boundary conditions (JTBl) and 
fTSOj) are fulfilled. Optimal third order polynomial protocols A*(r/t) are shown in Fig. [1] 
for different transition times t. A short transition time requires pronounced peaks with 
increasing height for decreasing transition times. These peaks serve to accelerate the 
Brownian particle in the first part of the transition and to deccelerate the particle in 
the second part of the transition. 

Another convenient choice is a linear protocol with superimposed delta peak 
singularities at beginning and end. The initial singularity allows for setting the initial 
mean squared momentum instantaneously to a value ^(0+) ^ z(0) = mT which can be 
tuned to meet the boundary condition on w(t) = w*(t). The second delta singularity at 
the final time t allows for setting w(t) =0. 

Note that, similar to the first case study, the minimal work (T2~T1) does not depend on 
the transition time t. Naively, in the limit t — > 0, one would expect the result obtained 
for an instantaneous jump of the protocol W iP = (A/ — A,)/2 which is substantially 
larger than the optimal work. However, the possibility to use a protocol with increasing 
absolute values of A (and thus increasing forces) for decreasing transition times t leads 
to the singular limit at t — > 0. For any fixed shape of the protocol X(r/t), the work 
approaches P0 P with t — > 0. 



Optimal protocols for Hamiltonian and Schrddinger dynamics 



7 



5 r 

4 - 

3 - 

2 j 

- 

-1 - 
-2 - 




-3 





0.2 



0.4 



0.6 



0.8 



1 



r/f 



Figure 1. Optimal third order polynomial protocols A*(r/t) for the harmonic 
potential with time-dependent stiffness, see equation (|12|) . with m = 1,T = l,k = 
l,Xt = 2,Aj = 1 for different transition times t. For short transition times, the 
optimal protocol shows a pronounced positive peak in the first part of the transition 
and a pronounced negative peak in the second part of the process, thereby regaining 
the invested work. 



2.3. Anharmonic quartic potential 

In both harmonic case studies, the optimal work W* is given by the adiabatic 
work W /ad , independent of the total transition time t. The optimal protocols A*(r) 
are highly degenerate. It is interesting to see whether these (unexpected) features 
persist for anharmonic potentials where the optimality conditions cannot be cast into 
only two conditions for the moments of the probability distribution. However, it is 
computationally quite difficult to determine optimal protocols in the anharmonic case. 
For a case study, we use the quartic potential 



with T = 1, Aj = 1, and Xf = 2 together with a Fourier ansatz for the optimal protocol 



where we omit higher odd frequencies because their symmetry is not suitable to produce 
the peak structure found in the two harmonic cases. We next optimize the parameters 
ak for a minimal mean work which we calculate as the average of work values from 
25000 randomly sampled initial values for xq and Pq. For the minimization, we use a 
standard Mathematica algorithm. Note that the optimal work values obtained in this 
numerical procedure are upper bounds for the true optimal work since we have used a 
finite number of parameters. 



V(x,X) = ]xx i 



(22) 



A*(r) = 1 H h a sm(irr/t) + V] a*, sin(2fc7rr/t) 

* k=i 



(23) 
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The adiabatic work can be calculated using the microcanonical distribution for 
given energy E and given A. The change of internal energy then is given by 

dX \ dX /micro ZJ pm (E-V(X,X)) dX 

with the microcanonical partition function 

dx . 1 (25) 

<j2m(E-V(x,\)) 

The evolution equation ([21]) is an ordinary differential equation for E(X). For the 
potential V(x, X) = Ax 4 /4, we obtain 

which has the solution E = EqX 1 ^ 3 with initial energy Eq. Averaging over initial 
conditions a;o,Po according to the appropriate initial Boltzmann distribution ([3]) then 
yields 

W = 3 ( ^-™ g 0.1949. (27) 

In Fig. we compare optimal work values with the adiabatic work. Even though 
it looks like convergence may be reached already by the small number of Fourier modes, 
we cannot decide this question conclusively. Particularly for very small transition times, 
optimal protocols feature extreme protocol values which may require a large number of 
Fourier modes. Moreover, we cannot exclude the possibility that the minimum search is 
trapped in a local minimum preventing the convergence towards the global minimum. 
Therefore, we cannot exclude the possibility that the adiabatic work can be reached in 
a finite time also in this case study. 

More likely, however, is the following scenario: The optimal work is discontinuous 
at t — 0. For any transition time, work values well below the work for an instantaneous 
jump Wi p = 0.25 but above the adiabatic work can be obtained. Clearly, even for very 
fast transitions compared to the characteristic oscillation time t c h ar — 3, the optimal 
work is well below the work for a linear protocol which has the limit W iV = 0.25 for 
t — > 0. The optimal protocol again shows pronounced peaks for short transition times t, 
see Fig. [2b. Since we use an ansatz with a finite number of parameters, it is impossible 
to decide whether the optimal protocol is degenerate. 

Even though we could not decide whether the adiabatic work can be reached in a 
finite time, we can show that free energy differences as determined via the Jarzynski 
relation indeed improve for an optimal protocol compared to a linear protocol, see Fig. 
[3]and its caption for details. 



3. Quantum Systems 



In this section, we investigate optimal protocols for quantum system. Fortunately, 
for harmonic systems which are comparable to the first two case studies from the 
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Figure 2. Optimization results for the anharmonic potential ([22| . (a) Optimal 
work from the Fourier ansatz (|23|) as a function of the transition time t for different 
number of Fourier modes n compared to the work W hn obtained for the linear protocol 
A lln = l+r/t and the adiabatic work W ad ~ 0.1949. (b) Optimal protocols as obtained 
from the Fourier ansatz (12311 with n = 3 for different transition times t. 
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Figure 3. Comparison of free energy estimates for the anharmonic potential ([22]) for (i) 
a linear protocol and (ii) the optimized Fourier protocol (n = 3) for t = 0.5. The data 
were obtained from Langevin simulation of 2 ■ 10 6 trajectories for each protocol with 
m = 1, T = 1, and A, = 1, A/ = 2. (a) Distribution P(W) of work values W for the two 
protocols. The columns show the free energy difference AF ~ 0.1733 and the mean 
work values (W«) ~ 0.2344(±0.0003), (WW) ~ 0.2025(±0.0002). (b) Histogram 

for the free 



(l//3)ln J2? =1 eM-m)/N 



of 2 • 10 5 Jarzynski estimates AF est ; 
energy difference obtained from N — 10 single trajectory work values W-j each. The 
mean squared error (MSE) of these estimates consists of two parts: the systematic error 
(bias) B = (AF cst ) — AF and the statistical error a — ^/Var(AF cst ). The columns 
show the free energy difference and the mean value of the estimates obtained from 
the two protocols. Since the bias (fiW = 0.0038(±0.0002), BM = 0.0018(±0.0002)) 
can be neglected for both protocols, the MSE is dominated by the statistical error 
(cr« = 0.0843, er(") = 0.0665) which is larger by 27% for protocol (i). 
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previous sections, work distributions p(W) have recently been calculated for Schrddinger 
dynamics [22j [23] . By using their results, it is easy to derive optimal protocols since in 
both cases, the work distributions (and therefore also the mean work) depend on a single 
"rapidity" parameter which vanishes for quasistatic transitions. In the last subsection, 
we study the two level system of a single spin in a time dependent magnetic field. If we 
impose that the absolute value of the magnetic field stays fixed, the adiabatic work can 
only be reached in a finite (non-zero) time. 



3.1. Moving harmonic potential 

We consider the quantum analogue of a moving harmonic potential 

2 

mil) 

l/(x,r) = T (x-A(r)) 2 (28) 
with Hamiltonian 



H = huj 



muo 2 . , 2 



+ -^m 2 (29) 



with creation and annihilation operators and a, Planck constant h, frequency ui and 
a time dependent minimum of the harmonic potential A(r) which is changed from an 
initial value A(0) = to a final value X(t) = A/ during a given finite time t. It has 
recently been noted that the mean work is given exactly by the classical expression for 
this case study (2?j. Thus, the optimal protocols as discussed in Sect. 12.11 also apply 
to the quantum case. However, the work distributions are different for classical and 
quantum systems. It has been shown (23J, [28] that the work distribution p{W) for an 
initially thermalized ensemble depends only on the absolute value \z\ of a "rapidity" 
parameter 



Wiri^M^ < 30) 

which vanishes for adiabatic transitions with t — > oo. However, even for transitions in a 
finite time t, there exist appropriate protocols A(r) which lead to a vanishing "rapidity 
parameter" z and thus to a vanishing mean work W* = W = 0. In addition to the 
optimal protocols from Sect. 12.11 it also suffices to consider the simple trigonometric 
protocol 

(An 2 n 2 -t 2 cu 2 )sin( 2n7rr ^ 



T 



t j 



t 2mrt 2 uJ 2 



(31) 



which can be easily shown to yield z = for any integer n. Again, the optimal protocol 
is highly degenerate. This degeneracy arises because the continuous protocol A(r) with 
an infinite number of parameters has to satisfy only two boundary conditions A(0) = 0, 
X(t) = Xf and the adiabaticity condition z — 0. 
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3.2. Harmonic potential with time dependent frequency 

A similar argument can be applied to the quantum analogue of the harmonic potential 
with time-dependent stiffness mA(r) with Hamiltonian 

= p^ ™ X (t)x 2 . (32) 

The control parameter A is changed time-dependently from an initial value A(0) = A, to 
a final value A(t) = Xf during a finite transition time t, which corresponds to a time- 



dependent oscillator frequency oo(r) = yA(r). It has been shown |22j that the work 
distribution and therefore also the mean work only depends on a single parameter 



Q = —^[\ % {x f X(t) 2 + X(tf 



+ {\ f Y(t) 2 + Y(t) 2 )\ -1 (33) 
where X(t) and Y(t) are solutions for the classical harmonic oscillator equation 

X(t) + A(r)X(r) = (34) 
with initial conditions 

X(0) = , X(0) = 1 , F(0) = 1 , Y(0) = 0. (35) 

The dynamics of X(r) and Y(r) is constrained by their time evolution with common 
protocol A(r) which yields the relation 

XY = YX. (36) 

Addition of the term XY on both sides and subsequent integration leads to 

XY = XY + 1 (37) 

where the integration constant is determined by the initial conditions fl35|) . Using 
standard techniques, we perform a minimization of Q under the constraint f|3T|) with 
respect to X(t), X(t), Y(t), and Y(t). The optimality conditions then are given by 



X*(t) = —=L==J y / A^7- XfY*(t) 2 

X*(t) = \I^Y*(t) (38) 
V Aj 

where Y*(t) can be chosen arbitrarily. The minimal value of the parameter Q then is 
Q* = and the optimal work is given by the adiabatic work which can be calculated 

as 

W* = W ad = -h(u f - Wi) coth(/3^/2) (39) 



where Uf = yA/ and Ui = Again, only a finite number of optimality conditions 

have to be fulfilled while the continuous protocol A(r) corresponds to an infinite number 
of parameters. The optimal protocol A*(r) thus is highly degenerate and the adiabatic 
work can be reached in any given (arbitrarily short) transition time. 
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3.3. Two-level system 

We next study a simple two level system consisting of a single spin in a time dependently 
rotating magnetic field B(r) = Be(r) where e(r) is a unit vector which is to be changed 
from an initial value e(0) = e z to a final value e(t) = —e z in a given finite transition 
time t, keeping the absolute value of the magnetic field B fixed. Here denotes the 
unit vector in i-direction, % e {x, The Hamiltonian of the system is given by 

H = (40) 

with Pauli matrices a and gyromagnetic ratio 7. In the basis of eigenvectors of a z , the 
initial canonical density matrix is given by 

p( ) = I ( e ^/ 2 I +) ( + | + e -^ B ' 2 \-) (-1), (41) 

where | +) is eigenvector of cr z with eigenvalue +1, | — ) is eigenvector of a z with 
eigenvalue —1 and Z is the canonical partition function. The work in such a finite 
time transition then is given by 

W = (E) t - (E} = tr{H(t)p(t)} - tr{tf(0)p(0)} (42) 

where the time evolution of the density matrix p(r) is given by the Liouville-von- 
Neumann equation 

d T p{r)= l -[p{r\H{r)]. (43) 

Quite generally, due to the adiabatic theorem in quantum mechanics, the 
occupation probabilities of the states do not change for quasistatic driving. The density 
matrix at the end of an adiabatic transition is therefore given by 

Pi = ^(e-^ B/2 |+) ( + | + e^ B ' 2 \-) (-1). (44) 

Thus, initial and final energy are equal tr{iJ(0)p(0)} = tr{H(t c )p(t c )} and the adiabatic 
work is zero W ad = 0. 

We next try to identify optimal protocols for such a transition. Generally, any 
pure spin state corresponds to a Bloch vector of unit length. Applying a magnetic 
field B leads to a precession of the Bloch vector around the axis of the magnetic field 
with angular frequency u L = jB. Therefore, after an appropriate time t c = h/ul, the 
direction of a Bloch vector perpendicular to the magnetic field is reversed. During the 
adiabatic evolution from initial state p to the final state pi, the initial Bloch vectors in 
e z and — e z directions are also reversed. In order to reach the adiabatic work in a total 
transition time t = t c , it therefore suffices to choose the protocol 

!e z for r = 
e x for < t < t c (45) 
-e z for t = t c 

which corresponds to an initial jump to a (constant) magnetic field perpendicular to the 
^-direction, i.e. B = Be x , and a final jump to the magnetic field B = —Be z . At the 
end of the transition, the density matrix is given by 

P(tc) = Pi (46) 
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and the mean work thus is equal to the adiabatic work W* = W /ad = 0. The adiabatic 
work thus can be reached in the finite time t c . In order to reach the adiabatic work 
for any larger transition time t > t c , it suffices to choose the protocol e*(r) for times 
t < t c , keeping the magnetic field B(r > t c ) = —Be z unchanged for all times r > t c . For 
transition times t < t c , the adiabatic work cannot be reached since there is no possibility 
to obtain the density matrix pi from the given dynamics in a time t < t c . 

However, if we had not constrained the magnetic field to have a given absolute value 
B but allowed for any value of the magnetic field during the transition, we would have 
found that the adiabatic work can be reached in any transition time t. This is obvious 
because the spin precession frequency is proportional to B and thus any precession 
frequency can be achieved in such a setup. For a given transition time t, the magnetic 
field then can be tuned to yield a precession frequency ojl = n/t which leads to W* = 0. 

More generally, if one allows any form of the Hamiltonian, it is obvious that a 
work value arbitrarily close to the adiabatic work can be reached in any given (short) 
time t for the following reason. Assume that the work W\ can be reached in a 
slow transition with dynamics ^\{t) for the time dependent Hamiltonian H\{t) with 
t\ ^> t. The Hamiltonian H(r) = (ti/t)Hi(t\r/t) then leads to an equivalent dynamics 
\P(t) = ^>\{t\r/t) and therefore to an equivalent final mean internal energy at time t. 
Thus, the work W\ (which is arbitrarily close to W ad for long times ti) can be reached 
in the (arbitrarily short) time t. Note, however, that this argument is no longer true if 
one allows only for any form of the potential V(x, r) since the kinetic energy (and thus 
the dependence on the momentum p) in the Hamiltonian H(t) is given and thus cannot 
be rescaled in this manner. 

4. Conclusion and Perspectives 

We have investigated optimal protocols for a minimal mean work for both Hamiltonian 
and Schrddinger dynamics with canonical initial conditions. This extends our previous 
studies of optimal protocols for overdamped and underdamped Langevin dynamics 
[241 [25] where we found unexpected jumps in the overdamped case and even delta 
peaks in the underdamped case. 

For Hamiltonian dynamics and harmonic potentials, the present optimization shows 
two unexpected features which do not occur for Langevin dynamics: (i) Independently 
of the transition time t, the optimal work W* is given by the adiabatic work obtained 
for an infinitely slow transition, and (ii) the optimal protocol A*(r) is highly degenerate 
since only two conditions on the moments of the probability distribution have to be 
fulfilled. 

For Hamiltonian dynamics in an anharmonic potential, we performed numeric 
calculations of the optimal protocol in a restricted state space. We found that the 
optimal work is well below the work for an instantaneous jump even for very fast 
transitions and we thus conjecture that the optimal work is quite generally discontinuous 
at t — 0. However, we could not decide from our numerical data whether the adiabatic 



Optimal protocols for Hamiltonian and Schrddinger dynamics 



14 



work can be obtained in an arbitrarily short time t also for anharmonic potentials. In 
order to decide this question conclusively, extensive numerics is necessary which will be 
left for future research. 

For both harmonic and anharmonic potentials, values very close to the adiabatic 
work can be reached with transition times well below the typical oscillation times. These 
results suggest that the efficiency of free energy calculations for Hamiltonian dynamics 
via the Jarzynski relation can be significantly improved by using the optimal protocol. 
This insight is even more relevant since thermodynamic integration leads to a biased 
estimate because of the inequality of adiabatic work and free energy difference. 

Quite generally, the question whether the adiabatic work can be achieved in an 
arbitrarily short time will depend on the dependence of the potential V(x, A) on the 
protocol A(r). In the harmonic case studies, any constraint on the maximal possible 
absolute value of A(r) for Hamiltonian dynamics would lead to optimal work values well 
above the adiabatic work for short transition times. Such a constraint could, e. g., be 
implemented by replacing A(r) = atanh(A(r)), which yields values between A m i n = —a 
and A max = a for any value of A. Consequently, the question arises whether the adiabatic 
work can be reached in an arbitrary short time using any type of intermediate potentials 
which does include not only those compatible with a given functional form V(x, A). It 
is obvious that, if one allowed also for a time- dependent mass m{r), the adiabatic work 
could be reached in any given time, in complete analogy to the discussion in the quantum 
case. 

For quantum systems in an initially thermalized state, minimal work protocols show 
the same features as in the classical case for harmonic potentials. For a two level system 
which could serve as a paradigm for a finite quantum system, it depends on the allowed 
parameter space for the "potential", i. e. the magnetic field, whether the adiabatic work 
can be reached in an arbitrarily short time. Generally, the adiabatic work can be reached 
in arbitrarily short time if any functional form of the time-dependent Hamiltonian H(t) 
is allowed. 

For most practical purposes, the range of the control parameter will be limited by 
the experimental conditions. For quantum systems, refined strategies for the optimal 
control of the time evolution of wave functions under such constraints have been 
developped [29j ESS EH [32], [33]. These techniques may be generalized in order to specify 
minimum work protocols for quantum nonequilibrium transitions in complex systems 
with additional constraints on the possible values of the control parameter A(r). 
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